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The low energy spectrum of a spin chain with OSp( 3|2) supergroup symmetry is studied 
based on the Betlie ansatz solution of the related vertex model. This model is a lattice 
realization of intersecting loops in two dimensions with loop fugacity z = 1 which provides a 
framework to study the critical properties of the unusual low temperature Goldstone phase 
of the O(N) sigma model for TV = 1 in the context of an integrable model. Our finite- 
size analysis provides strong evidence for the existence of continua of scaling dimensions, 
the lowest of them starting at the ground state. Based on our data we conjecture that 
the so-called watermelon correlation functions decay logarithmically with exponents related 
to the quadratic Casimir operator of OSp(3\2). The presence of a continuous spectrum is 
not affected by a change to the boundary conditions although the density of states in the 
continua appears to be modified. 


I. INTRODUCTION 

This paper is concerned with the study of the finite-size properties of a solvable two-dimensional 
vertex model based on the five-dimensional representation of the OSp(3\2) superalgebra. This 
system has a close relation with a particular Lorentz lattice gas used to model the diffusion of 
particles through randomly placed obstacles on the square lattice [1]. In this cellular automata the 
particle moves along the bonds of the lattice and is scattered according to scattering rules fixed a 
priori once it reaches a given node. Here the scatterers are constituted of mirrors tilted right and 
left, i.e. by ±^, with respect to the lattice When the particle collides with a mirror it will 

turn right or left, depending on the orientation of the latter. In the absence of a mirror at a node 
the particle passes the node on a straight path. The corresponding scattering rules are depicted 
in Figure [lj Amplitudes w%, W 2 and W 3 represent the fraction of right and left mirrors and node 
vacancies on the lattice, respectively. The kinetic properties of this lattice Lorentz gas have been 
investigated by numerical simulations where an anomalous diffusive behavior was observed flj. In 
the case of partially occupied lattice by mirrors (in 3 7 ^ 0 ) the fractal dimension of large trajectories 
was argued to be dj = 2 with the presence of logarithmic corrections mm- 
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FIG. 1. Scattering rules for right mirrors, left mirrors and empty sites. The different vertices appear with 
probabilities w 1 , W 2 , W 3 respectively. 


An alternative interpretation of the paths of particles in this lattice gas is in terms of the 
degrees of freedom of an intersecting loop model with periodic boundary conditions in both lattice 
directions imposed. In Figure [2] we show the three possible configurations of a node of such loop 
model together with their corresponding probabilities. Weighting every closed loop in a given state 



FIG. 2. The possible configurations of a node in the intersecting loop model and their associated Boltzmann 
weights w\, w 2 , W3. 


by a fugacity z the object of our study is the partition function, 

Z= £ (1.1) 

loop configurations 

where mi, m 2 and m 3 are the number of weights w \ 1 W 2 and W 3 and A f denotes the number of 
loops on a such statistical configuration. 

For the case of fugacity z = 1 this intersecting loop model can be reformulated as the OSp{ 3|2) 
supersymmetric vertex model which is the subject of this paper. The integrability of this model 
provides a framework for the study of its critical behaviour and - exploiting the equivalences listed 
above - of the peculiar properties observed in the Lorentz lattice gas. In fact, the finite-size analysis 
of the lowest excitation of the OSp{ 3| 2) superspin chain found the respective critical exponent to 
be very small |Tj . This was taken as an indication for the presence of a zero conformal dimension 
on the spectrum implying the superdiffusive behaviour ( df = 2) predicted for the Lorentz gas. 
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In the context of the loop model, it has been argued that this behaviour signals the existence 
of an unusual critical phase of intersecting loops: in two dimensions the crossing of loops, W 3 7 ^ 0 , 
is a relevant perturbation to the low temperature dense loops phase |Tj. As a consequence the 
long distance behaviour of correlation functions in the loop model with fugacity of closed loops 
z = N < 2 should be that of the Goldstone phase of the O(N) sigma model. For integer N this 
regime can be described in a supersymmetric formulation of the field theory in terms of m bosons 
and 2n symplectic fermions, N = m — 2 n. The OSp(3\2) vertex model is one from a class of 
integrable lattice regularizations of these models. Their central charges are c = N — 1 as expected 
in the Goldstone phase [Tj. 

There has been a series of attempts towards the identification of the other characteristic feature 
of this Goldstone phase starting from integrable lattice models, i.e. a finite density of vanishing 
critical exponents. The existence of a continuous spectrum of conformal weights has been estab¬ 
lished in several staggered superspin chains ishei. including a model based on the four-dimensional 
representations of U q {sl{2\\)\ alternating with their duals which, in the self-dual case, is isomorphic 
a deformation of the OSp( 212) chain relevant to the loop model with fugacity N = 0 {llj . Further 
details of the spectral properties of these models have been uncovered when it was realized that 
hidden within the zero charge sector of this superspin chain there exists a staggered six-vertex 
model which already displays a continuous low energy spectrum. For the latter strong evidence 
has been accumulated that the effective theory describing the low energy excitations of the model 
is the £L( 2 ,M)/t/(l) sigma model at a level related to the anisotropy [T2TiTT7j . We note, however, 
that the focus of these studies has been on the anisotropic deformations of the vertex models: while 
it has been established that the isotropic N = 0 models are on the boundary of the critical region, 
the question of the critical properties at the isotropic point itself has not been addressed. 

In this paper we want to return to the case of N = 1 as described by the supersymmet¬ 

ric vertex model. As argued above, the critical properties of this model are those of the proposed 
Goldstone phase of the 0(1) sigma model and can be studied based on its solutions by means of 
the algebraic Bethe ansatz. By means of an extensive finite-size study of the model we accumulate 
ample evidence for the existence of continua of critical exponents with lower edges at scaling di¬ 
mensions X = 0, 1, 2, 4, 6 ,... . With the exception of the the ground state of the superspin chain 
all of the states considered show strong logarithmic corrections to scaling governed by the flow of 
the model to weak coupling in the Goldstone phase. This is complemented by the observation that 
many excitations which have energies oc 1 /L for small system sizes but disappear from the low 
energy spectrum as the system size is increased. Both of these features of the spectrum require 
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very large system sizes to be studied for a reliable identification of the low energy effective theory: 
here we consider lattices with up to 4096 sites. 

Finally, since recent numerical studies of the N = 1 intersecting loop model have emphasized 
the importance of boundary conditions for the long distance behaviour of correlation functions 
USE] we consider both periodic boundary conditions to the superspin chain and twisted ones 
depending on the fermion number. Comparing the results we find that the amplitudes of the sub¬ 
leading (logarithmic in the system size) finite-size corrections do depend on the choice of boundary 
conditions. 


II. THE INTEGRABLE LATTICE MODEL 


The statistical configurations of the Lorentz lattice gas or intersecting loop model mentioned 
in the introduction have a one-to-one correspondence with the generators of the braid-monoid 
algebra. This fact has been elaborated previously in the work [1] but for sake of completeness 
we have summarized this equivalence in Appendix |A| This algebra can be used to built solvable 
models and it turns out that integrability is assured when the weights are parameterized as 


W\ = w 0 , w 2 = 


w 0 X 


w 3 = w 0 A 


( 2 . 1 ) 


1/2-A’ 

where A is a free spectral parameter and wq is an arbitrary normalization. This scale can be chosen 
to interpret the weights as probabilities but in our context we set it to unity. Note that all three 
weights are positive for 0 < A < 

Furthermore, it has been shown that this integrable manifold can be realized in terms of a stan¬ 
dard local vertex model. Its bond states are constituted of three bosonic and two fermionic degrees 
of freedom realized in terms of the five-dimensional representation of the OSp{ 3|2) superalgebra, 
see Appendix [B] The possible configurations on a vertex with their respective Boltzmann weights 
are encoded in the i?-matrix 


«oiW = E 

a,b= 1 


“l(-l Y aPbe ab ® e ba + W 2 a ^ a J e aJ ® Hd + W * e ab ® e , 

c,d= 1 


(?) 

ab 


( 2 . 2 ) 


(k) 

where e y ab are the 5x5 Weyl matrices acting either on the auxiliary space for k = 0 or on the 
quantum space associated to the sites of a chain of length L for k = 1, ■ ■ ■ , L. The symbol p a 
denotes the Grassmann parities distinguishing the bosonic (p a = 0) and fermionic (p a = 1 ) degrees 
of freedom. The 5 x 5-matrix a a b is the basic ingredient to built an explicit representation for 
the monoid operator and its expression depends much on the grading order basis, see for instance 
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[18j. Here we will consider two specific Grassmann orderings in which the two U{ 1) charges of the 
05^(312) algebra commuting with the operator Roj(X) are organized in a way which is suitable to 
perform the Bethe ansatz analysis. This turns out to be the fbbbf and bfbfb basis ordering and 
the corresponding forms for the matrix a are 


a 


(fbbbf) _ 


( 0 0 0 0 l\ 

0 0 0 1 0 

0 0 10 0 

0 10 0 0 

\ -1 0 0 0 0 


a 


(bfbfb) _ 


( 0 0 0 0 1 \ 

000-10 
0 0 10 0 

0 10 0 0 

\ 1 0 0 0 0 ) 


(2.3) 


For the vertex model on the square lattice with Lx L vertices and periodic boundary conditions 
for both bosonic and fermionic configurations in the horizontal direction we now construct the 
vertex model row-to-row transfer matrix. This operator is given as the supertrace over the auxiliary 
space of an ordered product of L matrices Roj(X ), 
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T(X) = [Rol(X)Rol-i(X) • • • i2oi(A)] e 


(2.4) 


a= 1 


As a consequence of integrability the transfer matrix commutes for different values of the spectral 
parameter, [T(A),T(/x)] = 0, and therefore generates a family of commuting operators. As usual, 
a Hamiltonian with local (nearest neighbour on the lattice) interactions is obtained by expanding 


(2.4) around the point A = 0 where the transfer matrix becomes proportional to the translation 


operator. The resulting expression for the integrable OSp{ 3|2) superspin Hamiltonian is 


3 = 1 


£<-u 

a, 6=1 


PaPb 


Si) J3+ 1 ) 


"ab 


"ba 


+ 2 E a ab® cd 


le ac } e bd 


( 3 + 1 ) 


(2.5) 


a,b,c,d=l 

where periodic boundary conditions for bosonic and fermionic degrees of freedom are assumed. 

In terms of the transfer matrix the partition function of the vertex model is given by a supertrace 
of the L th power of T(A), now taken on the 5^-dimensional quantum space. For system sizes up 
to 4 x 4 we find by direct computation 

5 l - T 2 


z = ^(~ 1 f 1+P2 ''' +PLT ^W = 


k= 1 


1 +A + 


A 

1/2 - A 


( 2 . 6 ) 


in agreement with the triviality of the partition sum (1.1) of the loop model with z = 1. Here 
p\. - ■ ■ ,p/ are the Grassmann parities of the degrees of freedom composing a given fc-state of 
the Hilbert space. The partition function of any statistical model is dominated by the largest 
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eigenvalue of the transfer matrix. In the present case the contribution of a single vertex to the 
partition function is the sum of the three weights w 1 +W 2 +W 3 . Therefore, as long as the Boltzmann 
weights are all non-negative (i.e. in the regime 0 < A < 1/2), we conclude that the largest eigenvalue 
of the transfer matrix A max (X) for the lattice with linear dimension L is 

L 


A-max(^) — 


1 + A + 


A 


(2.7) 


1/2 - A 

The fact that there are no subleading (in L) corrections is a consequence of the grading of the 
states together with the properties of the OSp( 3|2) representations appearing in the Hilbert space, 
see Appendix [Bj 


We note that as an immediate consequence of the expression (2.7) for the largest eigenvalue of 


the transfer matrix the ground state energy of the superspin chain (2.5) is 


En — —3 L 


( 2 . 8 ) 


without finite-size corrections. 


Eq. (2.6) implies that the partition function can be normalized to Z = 1 by rescaling of the 


local Boltzmann weights. Note that this does not necessarily mean that the low-lying excitations 
in the spectrum of the transfer matrix are trivial. In general, we can only infer that the critical 
properties are governed by a conformal field theory (CFT) with central charge c = 0. Since the 


Hamiltonian (2.5) is a non-Hermitian operator the continuum limit is not expected to be described 


by a unitary c = 0 conformal field theory. 

In order to study the finite-size properties of the low-lying spectrum of this quantum spin chain 
we turn to its Bethe ansatz solution. 

III. THE BETHE ANSATZ 


The diagonalization of the transfer matrix can be carried out within the algebraic Bethe ansatz 
framework. The essential tools have already been discussed before [18] and here we shall present 
only the main results. It turns out to be convenient to re-scale the spectral parameter by the 
imaginary unit to bring the resulting equations into a canonical form for performing numerical 


analysis. Denoting by A(A) the eigenvalues of the transfer matrix (2.4) we find that they can 
be parametrized in terms of two sets of rapidities {A^}, a = 1,2. The functional form of the 
eigenvalues as well as the algebraic Bethe equations satisfied by the rapidities depend on the 
grading chosen. 
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A. Bethe ansatz in fbbbf grading 

The Hilbert space of the superspin chain can be decomposed into the irreducible representations 
(irreps) (p\ q ) of OSp( 3|2) appearing in the tensor product (0; of local spins, see Appendix [b| 
In the grading fbbbf a highest weight state of the (8 L — 4)-dimensional representation (L — 1; 
is used as pseudo vacuum (or reference state) for the algebraic Bethe ansatz. Thanks to the 
algebra inclusion 05p(3|2) D 51/(2) © 51/(2) all states can be characterized in terms of the two 
U{ 1) charges £ z , s z from the SU( 2) subalgebras. In the reference state they take values L and 0, 
respectively. Bethe states in a (p; q) multiplet are parametrized by (L — n\) complex rapidities A^A 
and (L — ni — 7 ^ 2 ) complex rapidities A:- , where 


ni=p+l, n 2 = 2q - 1 


(3.1) 


for all states except the OSp( 3|2) singlet (0;0) for which n± = n 2 = 0 (note that only 05p(3|2) 
irreps with integer p appear in the Hilbert space of the superspin chain). 

The corresponding transfer matrix eigenvalue is given by the following expression: 

A - xf + i/2 L ^ A - A$ 1} + i/2 L -" 2 A - A$ 2) - i 


a(a)=- [a /i - i] l n 


+[am l n 


n 


A - a[ 2) 


A X-Xf-i/2 f = \ X-Xf-i/2 

L _„ 2 _ a (2) _ •) _ A (2) + i/2) L-ni A _ A (l) _ i L-rn-n 2 A _ A (2) + - /2 

3 = 1 


n v ta ; + n 7 j hd n j 

(i) 


A 


(1/2-zA) 
(—i/2 + A) _ 


L L—n 1 


(A-A A) (A — A:- — i/2) 


i=i A -A 


i=i A - A f - i/2 


n 

3= 1 


A-A^-i 


A-A 


(i) 


The eigenspectrum of the Hamiltonian can be derived from this expression giving 


ehaW},© 2 ’})^- E 


,( 2 ) 


L—n\ 


3 =i [a5 1} ] 2 + 1/4 


(3.2) 


(3.3) 


The sets of variables {A-} are constrained by the Bethe ansatz equations which for the fbbbf 
grading are given by 

L .rn - n2 A (l) _ a (2 ) + . /2 

11 hi) ; 

k =1 


aJ 1} + i/2 
LA^-i/2 


A - AA - i/2 


3 = !,••• ,L-n u 


L _i:> Af - a) 1 ’ + i/2 Af - Af + i/2 


i H-Af-i/2 


n 


(3.4) 


xf-X^-i /2 


j = !,-■■ ,L-m-n 2 . 
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B. Bethe ansatz in bfbfb grading 


In this grading the Bethe ansatz uses a highest weight state of the (8 L — 4)-dimensional repre¬ 


sentation (0;L/2) as reference state. States in the sector (p; q) with (3.1) are now parameterized 


by {L — n 2 — 1) rapidities A^ and (L — ni — 77 - 2 ) rapidities In terms of these parameters the 

corresponding transfer matrix eigenvalue is 

L ~n 2 ~l X _ A (l) _ i/ 2 L-U2-1 x _ A (l) + i/ 2 L-^-n 2 A _ H2) + % 


A(\)=[\/i+i] L n 


=1 A — Ap 4- i/2 


-m L n 


j=i X — X^p+i/2 

L . L -2JlT n2 (A - xf - i) (A - A f + i/2) A - aJ 1} X - xf - 3i/2 ' 

+ [X,l] M (A-Af) (A-Af-i/2) f = , A-Af-i X-Xf-i/2 


n 


A- A 


( 2 ) 


A 


L L—n 2 —I 


n 


A - AV 


( 1 ) 


(3/2-HA) 

(i/2 - A) J A _ A (i) _ i 


(3.5) 


The expression for energies of the Hamiltonian in this grading differs from the previous one by an 
overall minus sign: 

L—U2 — 1 

,(‘h = — 1 

3 = 


L—U2 — 1 

£ ^ + 1/4 


(3.6) 


The Bethe equations for the rapidities A,-, a = 1,2 in the grading bfbfb are 

j = !,-■■ ,L-n 2 - 1, 


A« + i/2 n L 


LAf - i/2. 


L -" 2 A^-A^+i/2 


A ( 1 ) _ A (2 > — i 12 
k =1 l ! ^ 


i-«-! a ( 2) _ A (1) + i/2 _ (A f - A® - i/2) (A® - A® + i) 

M Af-A«-i/2 = II (A® - A® + WO®-A®-if J = 


,L - m — n 2 
(3.7) 


Note that in this solution the inclusion OSp(3|2) D OSp{ 112) becomes manifest: the second Bethe 
equations are exactly those of the OSp{ 112) invariant vertex model in the presence of inhomo¬ 
geneities M- 


C. Example: L = 2 


For L = 2 the spectrum of the superspin chain decomposes into the OSp{ 3|2) singlet (0; 0) and 


two twelve-dimensional multiplets (0; 1) and (1; |), see Eq. (B6). For these we can identify the 
corresponding Bethe root configurations in both gradings: 


1 We use the same notation for the Bethe ansatz rapidities for both gradings. Therefore, whenever specific root 
configurations are discussed, they need to be seen in the context of the underlying grading. 
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For the fbbbf grading the Bethe equations are (3.4), the corresponding energy is (3.3). The 
three energies are 

• (1; \) (the fbbbf reference state): 

L — n\ = 0 Bethe roots A^, L — n\ — n 2 = 0 Bethe roots A^ 2 ): the energy of this state is 
E = +2. 


( 0 ; 1 ): 


here we have n\ = 1, ri 2 = 1 according to (3.1). The only one finite solution to the Bethe 


equations (3.4) is A^ = 0 giving energy E = —2. 


(0;0): 

for the singlet ground state the quantum numbers are (rii , n-z) = (0,0), therefore there are 
two roots on each level, a = 1,2. The unique solution to the Bethe equations is A^ = 0 
(degenerate roots) and A*, 2 ] = ±z^. The resulting energy is E = —6, as expected from the 
general considerations in Section [TT| 


For the grading bfbfb the Bethe equations (3.7) have to be solved 


(0; 1) (the bfbfb reference state) : 

the number of Bethe roots is L — 77,2 — 1 = 0 and L — n\ — n -2 = 0, respectively. With this 
we find E = —2, as in the other grading. 

(!; i) : 

With ni = 2, m = 0 there is only one finite solutions A^ 1 ) = 0 giving E = +2. 

( 0 ; 0 ): 

the singlet ground state is again parametrized by two roots on each level, which have to 


satisfy (3.7) in this grading. It is straightforward to find the solution to be A^ = 

A) 2 = 0 (degenerate roots) giving the ground state energy E = —6. 

Note that the degeneration of Bethe roots in the ’flat’ ground state is a feature which has also 
been observed in other (super-)spin chains [HIE, 19] . 

IV. GROUND STATE AND LOWEST EXCITATIONS 


As discussed above the ground state energy of the superspin chain is exactly given by (2.8). 


For L = 2 the corresponding root configurations in the Bethe ansatze (3.4) and (3.7) have been 
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obtained in the previous section where they were found to be singular in the sense that Bethe roots 
may degenerate. For odd chain lengths L the situation turns out to be easier: here we find that 
the root configurations describing the (0; |) ground state are non-degenerate. In the grading fbbbf 
it is given by collections of (L — l)/2 pairs of complex conjugate rapidities 

^ER, (4.1) 

on each level a = 1,2. In the thermodynamic limit, L —> oo, the deviations from these ’strings’ 
become small. This allows to compute the ground energy density and the Fermi velocity vf of the 
gapless low lying excitations within the root density approach [2D] giving lini£_> 0 a Eo(L)/L = —3 
and vf = 2,ir, see also Ref. 1. 

As a consequence of conformal invariance the leading terms in the finite-size scaling of energy 
levels are predicted to be 1211122] 


7 TVjP 

Eo(L) — Le oo = ——— c + o 
o L 


En{L) - E 0 (L ) = 


2itvf 

L 


X n + O ( — 


(4.2) 


Here c is the central charge of the effective low energy theory. As a consequence of (2.8) we 
have c = 0 without corrections to scaling for the superspin chain. From the scaling dimensions 
X n = h + h appearing in the low energy spectrum together with the conformal spin, s = h — h, 
which can be read off the momentum of the corresponding state, we can obtain the conformal 
weights ( h , h) of the operators in this CFT. We note that the CFT with c = 0 is not unitary which 


may lead to subleading finite-size corrections vanishing as inverse powers of logT in (4.2). The 
identification of these terms is one of the goals of this work. 

From the exact diagonalization of the superspin chain Hamiltonian for small system sizes we 
find that among the lowest excitations of the superspin chain there are the minimum energy states 
in the symmetry sectors corresponding to the (0; q) irreps of OSp(3\2) with q > 1. The highest 
weight state in these multiplets reached in the first Bethe ansatz (grading fbbbf) is described by 


(L — m, L — TT-i — 7 ^ 2 ) = (L — 1, L — 2 q) roots of (3.4). It turns out, however, that for this class of 


states it is much easier to work within the bfbfb grading, Eqs. (3.7). Here the highest weight state 
has (L — ri 2 — 1, L — nl — n2) = (L — 2q, L — 2q) roots X'-'' arranged into collections of complex 
conjugate pairs, namely 


\(i) ( 1 ) , 37 


, (2) (2) , * 
4 


( a ) r- 


(4.3) 


Configurations of this type with 2 q even (odd) can appear for L even (odd). The numerical 
estimates for the scaling dimensions obtained from solutions to the Bethe equations are slowly 
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decreasing with the system size once L is sufficiently large. For the lowest level (q = 1) this 
has been observed already in in Ref. 00 where this was taken as evidence for the presence of a 
zero conformal weight in the spectrum. Later, this feature has been argued to be one of the 
characteristics of the low temperature Goldstone phase of the O(N) sigma-model with N < 2 
realized in the presence of loop crossings |7j: the latter act as a perturbation which breaks the 
symmetry from O(N), or OSp(m\2n) with N = m — 2n, to O Sp(m\2n) / O Sp(m — l|2n). In the 
present case with N = 1, and using the system size L as a long distance cutoff, the single coupling 
constant of the resulting sigma model on this supersphere is found to be 


9(7 ~ log(L/Lo) (4 ' 4) 

within a perturbative RG approach [23J [23] (note that log Lq has to be negative for g a > 0 as 
expected on physical grounds). With this as an input we extrapolate the finite-size data for the 
scaling dimensions assuming a rational dependence on 1/logL. Our results for the lowest (0;g) 
states are displayed in Figure |3j From these data we conclude that all dimensions X^ 0;g ) with finite 
q > 1 vanish in the thermodynamic limit, showing subleading scaling corrections proportional to 
1/ log L. For q > 2 the fact that X( 0;g ) —> 0 can be obscured by the latter for quite large L rendering 
a finite-size analysis based on small system sizes impossible. 

To determine the amplitude of these subleading terms we note that this class of (p = 0; (^-states 
can be extended to include the ground states of the superspin chain, i.e. the singlet (0; 0) for L 
even and the quintet (0, for L odd. Since there are no finite-size corrections to the energies of 
the ground states (with q = 0 and q = | for L even and odd, respectively) we conjecture that the 


amplitudes of the logarithms are related to the quadratic Casimir (Bl) of OSp( 3|2) as 

q(2 q-l) 


X(p ;q) — 0 + 


1 3 


log(L/L 0 )’ * 2 ’ 2 ’ (4 ' 5) 

Comparing this conjecture with our numerical data we find very good agreement, see Figure [3](b). 
The lattice sizes L considered here do not allow for a readable estimate of the non-universal scale 


logLo though. In the context of the loop model the amplitudes (4.5) determine the long distance 
asymptotics of the ’watermelon’ correlation functions Gk(r) measuring the probability of k loop 
segments connecting two points at distance r : these correlators - two-point functions of the so- 
called /c-leg operators - vanish with a power of 1/logr, i.e. 


Gk(r) ~ 1/ (logr) Qfe . 


(4.6) 


The amplitudes conjectured for the periodic superspin chain imply that the exponents aj~ are to 
be taken from the set {2q{2q — 1), q = 1, 3/2, 2,...}, i.e. the eigenvalues of the quadratic Casimir 
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l/log(L) 


FIG. 3. Scaling dimensions X ( 0;9 ) extracted from the finite-size behaviour of the energies of the lowest states 
in the sector (0 ;q). These states have conformal spin 0. Data obtained for the superspin chain with even 
(odd) number L of sites are represented by filled (open) symbols. The dashed lines show the results of an 
extrapolation of the finite-size data (for system sizes up to L = 4096) assuming a rational dependence on 


1/log L 1 red circles indicate the results obtained within the string approximation (4.7) -(4.91. In (b) the 


data for the scaling dimensions are shown with the X ( 0;g ) rescaled with the conjectured amplitude (4.5). 


(Bl) of OSp( 3|2). Identifying the (0; q ) primary with the k = 2q-leg operator this agrees with RG 
calculations and numerical results for G^r) and G 4 (r) [17]. As will be seen below, however, the 
asymptotic behaviour of the watermelon correlators is likely to depend on the boundary conditions. 


We further note that the reference state in the second Bethe ansatz is a highest weight state in 
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the unique (0; L/ 2) multiplet of the OSp{ 3|2) superspin chain of length L. This indicates that the 
spectrum of the (0; q) states considered here extends from the ground state energy Eq = —3 L to the 


energy E( 0; L/ 2) = — L. This suggests that the critical dimensions (4.5) actually form a continuum 
starting at A^o) = 0. Additional support for this suggestion comes from the observation that for 
q sufficiently large the imaginary parts of the Bethe roots in the bfbfb grading are exponentially 


(in L ) close to the hypothesis (4.3). In addition the real parts of the (L — 2q)/2 pairs on each of 


the two levels take values very close to each other, i.e. p 


(i) 




( 2 ) 


In this situation the second 


set of the Bethe equations (3.7) is automatically satisfied while the first level ones become, after 
taking their logarithm, 

(L-2q)/2 

L [^5/4(Mj) - Vh/4 (h)\ = -ZnQj + Y [V>3/2(Mj - «fc) + V’l(Mj - «fc) - V>l/2 {Vj ~ Mfc)] ■ 

k+j 

(4.7) 

Here 3 p a (x) = 2 arctan (x/a) and the numbers Qj define the many possible branches of the logarithm 
being given by the expression, 


1 


1 


1 


Qj = -t(£ - 2 g) + j - -, j = -(L - 2 q) . 


(4.8) 


4 v J 2 2 

Within this approach the eigenenergies corresponding to this state can be obtained using the 
following expression 


(L—2q)/2 

E({ H }) = -L+ 2 Y 

3 = 1 


1 


L(Ml) 2 + (5/4) 2 (Ml) 2 + (1/4) 2 J ' 


(4.9) 


Comparing the numerical solution of these ’string’ equations with the those of (3.7) we conclude 


that the finite-size energies in the sectors (0; q) with q > 2 are reproduced by Eqs. (4.7)—(4.9), see 
Figure [3](a). In this formulation the root density approach can be applied to compute the finite- 
size energies. In this approach we find again that for L —> oo the conformal dimensions are zero, 
independent of q. This observation provides an additional analytical support to the existence of a 
continuum of zero conformal weights starting at zero. 


V. OTHER EXCITATIONS 

As for the sectors (0; q) in the previous section we have identified the Bethe configurations 
corresponding to the lowest energy states in the sectors (p: q) with p ^ 0. In Figure [4] we present 
the corresponding scaling dimensions for p = 1. Again, the energies show strong logarithmic 
corrections to scaling which are dealt with in the extrapolation by assuming a rational dependence 












14 



_,_ l _,_l_ i _l_ 

0 0.2 0.4 0.6 

1/1 ogL 


FIG. 4. Same as Fig. [3^ a) but now for the lowest states in the sector (l;q). These states have conformal 
spin 1. 

of the finite-size data on 1 / log L. For all states which we have considered the scaling dimension 
extrapolates to = 1. All of these states are found to carry momentum ±2n/L, from which 

we conclude that the corresponding fields have spin 1. 

The lowest energy levels in some sectors (p; q ) with p > 2 are studied in an analogeous way, 
see Figures [5j [6j and [7j The extrapolation gives scaling dimensions X( 2]q ) = 2, X ( 3;? ) = 4, and 



FIG. 5. Same as Fig. [3^ a) but now for the lowest states in the sector (2 ;q). These states have conformal 
spin 0. 
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FIG. 6. Same as Fig. [3^ a) but now for the lowest states in the sector (3 ;q). These states have conformal 
spin 2. 



FIG. 7. Same as Fig. [3^ a) but now for the lowest states in the sector (4;g). These states have conformal 
spin 0. 


= 6 as lower edges of continua of states with different q. The low lying states in sectors 
with p = 2,4 have zero momentum while the states in the sectors (3; q ) have spin 2. Note that the 
presence of these levels in the spectrum of the superspin chain of length L, requires the selection 
rule L — p — 2q E 2N to be satisfied. 

In Table [T] we have collected the discrete parts of the conformal weights of primary operators 
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in the CFT as identified from our numerical studies. Based on these data we conjecture that the 


TABLE I. Discrete part of the scaling dimensions X , conformal spin s and corresponding conformal weights 
(h, h ) obtained from the lowest states observed in the symmetry sectors (p; q) of the superspin chain. 


sector (p; q) 

X = h + h 

s = h — h 

(h,h) 

(0 ;<?) 

0 

0 

(0,0) 

(!;<?) 

1 

±1 

(L0), (0,1) 

M 

2 

0 

(LI) 

(3; q) 

4 

±2 

(3,1), (1,3) 

(4; q) 

6 

0 

(3,3) 


spectrum of conformal weights is given by hk = \k(k + 1), k = 0,1, 2 ,... and the lowest levels in 
the (p; q) sector of the spectrum of the superspin chain correspond to operators with conformal 
weights 


(h,h) = 


i h P / 2 : hp/'i) for p even, 

(h(p± i)/2, h(p-_ fi)/ 2 ) for p odd. 


(5.1) 


Having identified the primary operators responsible for the lowest energy states we proceed to 
studying the finite-size behaviour of the other states appearing in the sector (0; 1): from the numer¬ 


ical diagonalization of the Hamiltonian (2.5) for the superspin chain with L = 6 sites we observe, 


apart from the (0; 1) primary with vanishing scaling dimension (4.5), two low lying excitations 


with spin ±1 and three low lying excitations with zero spin. We have identified the configurations 
of Bethe roots corresponding to all of these except one of the spin 1 states. Solving the Bethe 
equations for these configurations for larger system sizes and computing the corresponding scaling 
dimensions we find that they extrapolate to X = 2 for two of the spin 0 states, indicating a de¬ 
scended field of the (0; 1) primary operator. The scaling dimensions of the third spin 0 state and 
the spin 1 state, however, do not converge to a finite value under extrapolation assuming a rational 
dependence on 1/logL but rather disappear from the low energy spectrum, see Figure |8j 

This behaviour is not limited to excitations in the (0; 1) sector but we have found such levels in 
many other sectors, too. Let us note that the removal of low energy states present in small systems 
in the scaling limit has also been observed in one phase of the U q [sl{ 2| 1)] staggered superspin chain 


m ■ It is a direct consequence of the vanishing of a coupling constant, such as in (|4.4|) for the 


OSp(3\2 ) model, and therefore expected to be a generic property of realizations for field theories 
with a continuous spectrum of critical exponents as lattice models with a compact quantum space. 
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FIG. 8. Scaling dimensions X extracted from the finite-size behaviour of the energies of some low lying 
states in the (0; 1) sector together with their extrapolation assuming a rational dependence on 1/logL. The 
levels with X ~ 1 (1.75) for small systems (largest 1/logL) have conformal spin 1 and 0, respectively. 


From conformal field theory one expects a spin 1 descendent of the (0; 1) primary field with 
scaling dimension X = 1. This might be the other low energy spin 1 level observed in the L = 6 


spectrum. Since we have not able to identify a solution to the Bethe equations (3.4) or (3.7) 
reproducing the numerical value for the energy, we cannot confirm this expectation. 


VI. BOUNDARY CONDITIONS 


Here we would like to point out that peculiar finite-size behavior reported so far is also present 
if we had modified the boundary conditions such that the five possible states on the bonds of the 
vertex model are considered to be bosonic degrees of freedom. In this situation the vertex model 


transfer matrix is given, instead of (2.4), as the standard trace over the auxiliary space of the 
following product of operators, 

5 


T(A) = £ [I OL RoL{X)ioL-iRoL-i{X) ■ ■ ■ loi^oiW], 


( 6 . 1 ) 


a =1 


where the operator Iqj = Y^ a b=i(~l) PaPbe aa x pl a y s the role of a graded identity matrix. We 
recall here that this type of transfer matrix has been considered before in the case a vertex model 
based on the OSp( 112) superalgebra [251]. 

From the point of view of the spin chain this modification corresponds to imposing anti-periodic 
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boundary for the fermionic degrees of freedom which splits the (p; q) multiplets of the superspin 
chain into two subsets (p; q) e and (p; q) Q containing states with even and odd number of fermions, 
respectively. In the Bethe equations this twist is reflected by additional phase factors depending 
on the sector label of the conserved charges. In the fbbbf grading only the first set of the Bethe 


equations (3.4) is affected by this change of boundary conditions giving 

^ L L—n\—ri 2 \ (1) __ \ (2) , • /n 

= (-ir +i n V V 7 , 4 i. 


+i/2 




(i) 


- i/2. 


^ - 4 2) 


L r ^ 1 A f - \£ ] + i/2 _ L -^- n2 AS 2) - A® + i/2 
II TM To) 77 ~~ 11 

kj=j 


-- \ (2) x(l) ■/« 

fc=i A i — A k 


k=1 

"J ''fc 

xf-Xf-i /2 


-i/2 

/ = !,■■• , 7 - ni - n 2 . 


( 6 . 2 ) 


The corresponding energies of the Hamiltonian with these boundary conditions are given by 


Eq. (3.3), as before. 


It is evident from this construction that the finite-size spectrum will be different from that of 
the superspin chain only for sectors where the number of fermions n i is even. To see the effect on 
the scaling dimensions we have to take into account that the Bethe state of highest weight in the 
even fermion sector of the (p; g)-multiplet is parametrized by 


ni = p, U 2 = 2 q 


(6.3) 


rapidities in (6.2) for p even rather than (3.1) which still holds for p odd 


The ground state of the spin chain with twist is the unique singlet (0; 0) as for the superspin 


chain before. The corresponding root configuration consists of L/2 pairs of rapidities (4.1). Unlike 


the situation in the superspin chain the ground state of the twisted model has a strong finite-size 
dependence on L, see Figure [9](a) : the additional phase in the Bethe equations (6.2) leads to a 
central charge c = 3 as noted before in Refs. Q] and [7, In addition there are subleading corrections to 
scaling which turn out to much stronger than the l/[logL] 3 behavior usually observed in isotropic 
spin chains. We stress that this is even in contrast to the finite-size behavior found for other spin 
chains invariant by superalgebra such as OSp( 112) and OSp(2|2). 

In addition we have solved Eqs. ( |6.2| ) for some low energy excitations (p; q) e with even fermion 
number n\. Their root configuration contains n 2 real roots in addition to pairs as in the 
ground state. Our results show that the low energy spectrum of the spin chain with anti-periodic 
boundary conditions for the fermionic states shows a similar behaviour similar to that of the 
superspin chain, see Figure [9](b): the lowest excitations have fermion number m = 0 with a 
continuous spectrum starting at scaling dimension X e = 0. As discussed above, the energies of 
states with odd fermion n\ number are parametrized by the same rapidities as the corresponding 
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1/logL 


FIG. 9. (a) Central charge c and (b) scaling dimensions X extracted from the finite-size behaviour of the 
energies of the ground state and some low lying states (p; q) e with even fermion number of the spin chain 
subject to anti-periodic boundary conditions for the fermionic degrees of freedom. Dot-dashed lines in (b) 
indicate scaling dimensions extracted from the corresponding levels ( p ; q) a with an odd number of fermions. 


(p; q ) multiplet of the superspin chain. Here the scaling dimensions, however, have to be computed 


from (4.2) relative to the new ground state which leads to a shift X Q = X( p;g) +c/12 as a consequence 


of the different effective central charge. As in the superspin chain the subleading corrections to 
the scaling dimensions vanish as 1/logL. The amplitudes of these terms display a q-dependence 
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which is clearly different from (4.5), even for the X a as a consequence of the logarithmic corrections 
to the central charge. Our data do not allow to quantify these amplitudes though: (much) larger 
system sizes would be needed for an estimate which is beyond the methods used in this work. 
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Appendix A: Correspondence between the braid-monoid and the superalgebra realization of 

the intersecting loop model 

The scattering rules of the Lorentz lattice gas and the configurations of the intersecting loop 
are algebraically realized in terms of the generators of a braid-monoid algebra first introduced by 
Brauer [26]. The braids Bj correspond to the right mirrors, the monoids E j represent the left 
mirrors while the identity I j describes the absence of a scatterer. These operators acts on the sites 
j = 1, • • • , L of chain of length L and satisfy a number of algebraic constraints: 

• Braid relations 

I I > j : iB.i B ;: |B;B, : |. B 2 = Ij, BjBfc = B fc Bj for \j - k\ > 2 (Al) 

• Monoid relations 

EjEj-tiEj = E j, E j = zFij, EjEfc = E^.Ej for |j — k\ > 2 (A2) 

where z is a free parameter representing the fugacity in the loop model version. 

• Mixed relations 


BjEj = EjBj = E j, EjBj-tiBj = Bj-|-iBjEj-|-i 




(A3) 


This algebra admits the construction of a one parameter family of integrable models by means 
of the approach known as Baxterization m ■ The respective i?-matrix is built out in terms of a 
weighted linear combination of the generators, 


Rj = 1V3 lj + W 2 E j + u>i Bj . (A4) 

It turns out that the matrix Rj fulfills the Yang-Baxter equation provided that the weights w\, 
W 2 and W 3 are sited on the projective quadric, 

2 - z 

(wi + W2)w3 --— W1VJ2 = 0 . (A5) 

This manifold can be parameterized with the help of one spectral parameter A as follows, 


wo A . 

wi = w 0 , w 2 = - -tv -r , ^3 = wq\ 

1 — z 2 — A 


(A6) 


where wq is an overall normalization. Note that this factor can be chosen such that the weights 
configurations are interpreted as probabilities, i.e. wq = (z — 2 + 2A)/((’: — 2)(1 + A) + 2A 2 ). 
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It has been observed [1J that for integer 2 this algebra has a realization in terms of a finite 
dimensional representation of the superalgebra OSp(m\2n ) provided that 


z = m — 2n. 


(A7) 


In this formulation the braid operator becomes the graded permutation between m bosonic and 
2 n fermionic degrees of freedom, 

m+2n 


Bj = £ (-i) p ^ 


® e ba +1) > 


(A8) 


a,b= 1 


where p a is the Grassmann parity of the a-th degree of freedom assuming values p a = 0 for bosons 
and p a = 1 for fermions. 

Similarly, the monoid operator can be written as, 


m+2n 


E i = £ a abaJ- {3) - 0+1) 

a,b,c,d=l 


e ac <» %d 


(A9) 


where the matrix elements a a b are 

(—l) 1_Pa for a = b 


®-ab * 


+1 for a < b with p a = Pb = 1 and ^£ 


= n 


k=a -\-1 


(A10) 


—1 for a > b with p a = Pb = 1 and = n 

k=b -\-1 


0 otherwise 

in terms of the gradings p a . 

We remark that for some specific grading ordering the matrix a can be transformed into a block 


anti-diagonal structure. For a generic example see Ref. [T8l and the main text, Eq. (2.3), for the 
case of OSp( 3|2). 


Appendix B: Finite dimensional representations of OSp( 3|2) 

As a consequence of the algebra inclusion OSp( 3|2) D SU( 2) © SU(2) one can use the labels 
of 577(2) © SU(2) irreps to classify the basis states of an OSp{ 3| 2) irrep, see [28]: except for 
the trivial representation (0; 0) the latter are characterized by two integer or half integer numbers 
(p > 0; q > |). The quadratic and quartic Casimir operators of OSp( 3|2) in terms of these numbers 
are 

h = {pip + 1) + 2q(l - 2q )), I 4 = £ 2 (3p{p + 1) + 2(q + l)(2q - 3)). (Bl) 
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For example, the five-dimensional fundamental representation (0; |) of OSp(3\2) carried by the 
local spins in the superspin chain decomposes as 


(0; |) = (1 <S> 0) © (0 <8) |) 


(B2) 


into (£ <S> ^-representations of SJJ{2) © SU(2). States with integer (half-integer) spin in the second 
factor are bosonic (fermionic). Similarly, the first few of the other 05p(3|2)-representations with 
integer p (i.e. the only ones relevant for the OSp{ 3| 2) superspin chain) can be decomposed as 

(0; 1) = (0 <g> 1) © (1 <g> 0) © (1 <g> i), 

1 3 

(0; q > 1) = (0 (8) q) © (1 (8) q - 1) © (1 ® q - -) ® (0 (8) q — -), 

1 1 

(p > 1; -) = (p - 1 ® 0) ® (p + 1 <g> 0) © (p <g> -), 

(1; 1) = (1 <g> 1) © (2 <g> 0) © (1 ® ^) © (2 ® ^), 


(p > 1; 1) = (p<g> 0) © (p< 


1 , 


(p <g> 1) © (p - 1 <8> 0) © (p + 1 <g> 0) 


(B3) 


. 1 N . 1. 

© (p - 1 ®-) © (p + 1 <8>-), 

For p > 1 and q > | with p ^ 2q — 1 the “typical” OS'p(3|2)-representation (p; q) has dimension 
4(2p + 1)(4 q — 1) and is an eightfold pattern in the 577(2) © 577(2) decomposition 


(p; q) = {p ® v) © (p ® q ~ g) © (p - 1 ® Q ~ \) © (p + 1 ® q - \) 

3 

®(p®q-l)®(p-l®q-l)®{p+l®q-l)®{p®q-~). 


(B4) 


For p = 2q — 1 some of the matrix elements of the OSp( 3|2) generators vanish which allows to 


decompose (B4) into two “atypical” representations with a fourfold 577(2) © 577(2) decomposition 


(p; q = ^(p + 1)) = (p® q) © (p © q - ^) © (p + 1 <8> q ~ ^) © (p + 1 © q ~ 1) • (B5) 


Note that the atypical representations together with the irreps (0; ^), (1; 1) and the trivial rep¬ 


resentation (0; 0) cannot be uniquely specified by the eigenvalues of the Casimir invariants (Bl): 
for all of them I 2 = I 4 = 0. Therefore, they can appear as parts of reducible but indecomposable 
representations of 05p(3|2) which are present in the spectrum of the superspin chain as will be 
discussed below. 

Based on these relations the Hilbert space of the isotropic OSp{ 3| 2) superspin chain for small 
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L can be decomposed into 0,Sp(3|2)-multiplets 


- ind 


(B 6 ) 


( 0 ; |)® 2 = ( 0 ; 0 )©( 0 ; 1 ) © ( 1 ; |), 

( 0 ; = ( 0 ; \) © ( 0 ; |) 0 ( 2 ; 1 ) © 2 * [ 2 * ( 0 ;|) 0 ( 1 ; 1 ) 

(0; l)® 4 = (0; 0) © 6 * (0; 1) © (0; 2) © 6 * (1; ^) © 3 * (1; 

© 3 * (2; 1) © (3; ^) © 2 * [2 * (0; 0) © (1; l)] ind . 

By [.. •lind we denote the combinations of atypical representations in reducible but indecomposable 


ones. The decompositions (B 6 ) have been checked numerically. Note that the trivial and the 
fundamental representation of OSp(3\2) appear exactly once in the tensor product for even and 
odd L, respectively]^] Furthermore, all other multiplets appearing in (B 6 ) have even dimension 
and contain the same number of bosonic and fermionic states. This leads to the trivial partition 


function (2.6) of the OSp( 3| 2 ) vertex model. 

Reference states to be used in the algebraic Bethe ansatz of the superspin chain with L > 1 


sites, each carrying the representation (B2), can be highest weight states of (either of) the ( 8 L —4)- 
dimensional multiplets (L — 1; or (0; \L). 


2 Other instances of these representations appear as parts of the indecomposables, though. 
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